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decreases  significantly  with  well  number  up  to  at  least  three,  and  that  this  decrease  can  occur 
without  sacrificing  much  of  the  current  density  or  speed  of  the  devices.  The  minimum  room- 
temperature  shot  noise  factor  predicted  by  this  work  is  0.24  in  a  triple-well  structure. 

This  project  has  also  paved  the  way  for  a  new  type  of  “Fermi-electronic”  device  -  the 
resonant-plasmon-mediated  switch.  In  this  type  of  device,  the  quantum  well  of  a  resonant 
tunneling  devices  would  be  replaced  by  a  semiconducting  matrix  filled  with  a  3D  distribution 
of  semimetallic  nanoparticles.  Each  nanoparticle  is  capable  of  storing  charge  and,  therefore, 
in  providing  electrostatic  feedback  to  free  carriers.  Because  the  nanoparticles  are  not  directly 
in  the  path  of  the  charge  carriers,  it  would  be  possible  to  get  the  noise  suppression  of 
resonant-tunnel  diodes  but  over  a  much  broader  range  of  voltage  and  at  much  higher  current 
densities.  This  type  of  mechanism  should  be  capable  of  smoothing  out  thermal,  or  Johnson- 
Nyquist  noise,  when  the  device  is  biased  out  of  thermodynamic  equilibrium. 

As  a  first  test  of  resonant  plasmon-mediated  structures,  near-inffared  attenuation 
measurements  were  made  between  X  =  0.5  pm  (hv  =  2.48  eV)  and  10.0  pm  (hv  =  0.124 
eV)  on  three  GaAs  epitaxial  layers  containing  thin  ErAs  layers.  At  wavelengths  in  the 
vicinity  of  2.5  pm  well  below  the  GaAs  band  edge,  an  asymmetric  res|>nant  attenuation 
signature  was  observed  that  increases  in  strength  monotonically  with  the  ErAs 
nanoparticle  density.  The  peak  attenuation  coefficient  was  found  to  be  1 300, 1 850,  and 
2380  cm’’  for  ErAs  volumetric  flections  of  0.56,  0.85,  and  1.13  %,  respectively.  The 
attenuation  signature  is  fit  well  by  a  Maxwell-Gamett  composite-dielectric  model  that 
treats  the  ErAs  as  semimetallic  spherical  particles  having  carrier  concentrations  and  the 
corresponding  effective  masses  in  accordance  with  separate  band-structure  computations. 

The  resulting  absorption  peak  corresponds  to  the  fundamental-mode  surface-plasmon 
(i.e.,  Frohlich)  resonance  of  the  ErAs  nanoparticles  embedded  in  GaAs. 
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Abstract 


This  project  has  developed  a  formalism  for  computing  the  shot  noise  in  resonant¬ 
tunneling  structures  having  an  arbitrary  number  of  quantum  wells.  It  is  based  on  the 
electrostatic  feedback  effect  in  quantum-well  structures  and  is  the  first  known  shot  noise 
formalism  that  treats  the  electrostatic  and  quantum  effects  in  a  self-consistent  fashion. 
Simulation  results  are  analyzed  and  compared  to  experimental  results  for  an  InGaAs/AlAs 
double-barrier  structure,  which  is  the  fastest  resonant  tunneling  diode  ever  characterized  for 
its  noise  properties.  Shot  noise  reduction  and  enhancements  are  predicted  in  the  positive  and 
negative  differential  resistance  region,  respectively,  in  satisfactory  agreement  with 
experiment.  The  shot  noise  reduction  effect  is  then  studied  for  multiple-well  resonant 
tunneling  diodes.  It  is  found  that  in  properly  designed  structures  the  shot  noise  factor 
decreases  significantly  with  well  number  up  to  at  least  three,  and  that  this  decrease  can  occur 
without  sacrificing  much  of  the  current  density  or  speed  of  the  devices.  The  minimum  room 
temperature  shot  noise  factor  predicted  by  this  work  is  0.24  in  a  triple-well  structure. 

This  project  has  also  paved  the  way  for  a  new  type  of  “Fermi-electronic”  device  -  the 
resonant-plasmon-mediated  switch.  In  this  type  of  device,  the  quantum  well  of  a  resonant 
tunneling  devices  would  be  replaced  by  a  semiconducting  matrix  filled  with  a  3D  distribution 
of  semimetallic  nanoparticles.  Each  nanoparticle  is  capable  of  storing  charge  and,  therefore, 
in  providing  electrostatic  feedback  to  free  carriers.  Because  the  nanoparticles  are  not  directly 
in  the  path  of  the  charge  carriers,  it  would  be  possible  to  get  the  noise  suppression  of 
resonant-tunnel  diodes  but  over  a  much  broader  range  of  voltage  and  at  much  higher  current 
densities.  This  type  of  mechanism  should  be  capable  of  smoothing  out  thermal,  or  Johnson- 
Nyquist  noise,  when  the  device  is  biased  out  of  thermodynamic  equilibrium. 

As  a  first  test  of  resonant  plasmon-mediated  structures,  near-infrared  attenuation 
measurements  were  made  between  X  =  0.5  pm  (hv  =  2.48  eV)  and  10.0  pm  (hv  =  0.124 
eV)  on  three  GaAs  epitaxial  layers  containing  thin  ErAs  layers.  At  wavelengths  in  the 
vicinity  of  2.5  pm  well  below  the  GaAs  band-edge,  an  asymmetric  resonant  attenuation 
signature  was  observed  that  increases  in  strength  monotonically  with  the  ErAs 
nanoparticle  density.  The  peak  attenuation  coefficient  was  found  to  be  1300,  1850,  and 
2380  cm’’  for  ErAs  volumetric  fractions  of  0.56,  0.85,  and  1.13  %,  respectively.  The 
attenuation  signature  is  fit  well  by  a  Maxwell-Gamett  composite-dielectric  model  that 
treats  the  ErAs  as  semimetallic  spherical  particles  having  carrier  concentrations  and  the 
corresponding  effective  masses  in  accordance  with  separate  band-structure  computations. 

The  resulting  absorption  peak  corresponds  to  the  fundamental-mode  surface-plasmon 
(i.e.,  Frohlich)  resonance  of  the  ErAs  nanoparticles  embedded  in  GaAs. 

1.  introduction 

For  the  past  forty  years,  electronic  computers  have  grown  more  powerful  as  their 
basic  sub-unit,  the  transistor,  has  shrunk.  However,  the  laws  of  quantum  mechanics  plus 
the  limitations  of  materials  and  fabrication  techniques  soon  are  likely  to  inhibit  further 
improvement  in  the  performance  with  reduction  in  the  size  of  today's  bulk-effect 
semiconductor  transistors.  Researchers  have  projected  that  as  the  overall  size  of  the  bulk- 
effect  semiconductor  transistor  is  miniaturized  to  approximately  0.1  micron  and  below, 
the  devices  may  not  improve  significantly.  Thus,  in  order  to  continue  the  miniaturization 
of  integrated  circuits  well  into  the  next  century,  it  is  likely  that  present-day  micron-scale 
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device  designs  will  be  replaced  with  new  designs  or  new  device  types  that  take  advantage 
of  the  quantum  mechanical  effects  that  dominate  on  the  much  smaller,  nanometer  scale. 
Quantum-effect  nanoelectronic  devices  are  a  promising  alternative  for  continuing  to 
increase  the  density  and  speed  of  integrated  circuits. 

Among  the  numerous  nanoelectronic  approaches  proposed  (Fig.  1)  and 
demonstrated,  resonant  tunneling  are  a  strong  candidate  for  circuit  applications  due  to 
their  negative  differential  resistance  (NDR)  characteristics,  structural  simplicity,  relative 
ease  of  fabrication,  inherent  high  speed,  flexible  design,  and  versatile  circuit 
functionality.  There  are  also  good  practical  reasons  to  believe  that  the  resonant  tunneling 
diode  (RTD)  may  be  the  next  device  based  on  quantum-confined  heterostructures  to 
make  the  transition  from  the  world  of  research  into  practical  application.  Progress  in 
epitaxial  growth  has  improved  the  peak-to-valley  current  ratio  at  room  temperature  even 
beyond  that  required  for  many  circuit  applications.  This  is  what  distinguishes  the  RTD 
from  other  interesting  quantum  device  concepts  that  show  weak,  if  any,  desired 
phenomena  at  room  temperature. 

Moreover  the  study  of  resonant  tunneling  diodes  provides  valuable  knowledge  of 
the  quantum  aspects  of  electron  transport  in  mesoscopic  systems.  Tunneling  current  is 
becoming  a  more  interesting  issue  in  today’s  transistors.  Quantization  of  carriers  in  one, 
two  or  three  dimensions  is  a  promising  basis  of  device  design  in  the  future.  The  focus  of 
the  present  work  is  the  noise  characteristics  of  resonant  tunneling  diodes. 

As  semiconductor  devices  have  shrunk  in  size,  the  intrinsic  noise  is  more  than 
ever  a  limiting  factor  of  their  capabilities,  especially  in  nanoelectronic  devices.  Given  a 
specific  device  type,  the  traditional  approach  to  noise  control  and  suppression  has 
occurred  at  the  circuit  level  where  a  large  number  of  techniques  (e.g.  negative  feedback) 
have  been  developed  to  achieve  low-noise  circuits  with  garden-variety  devices.  The 
present  work  is  aimed  at  engineering  noise  suppression  in  the  semiconductor  device 
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Fig.  2  Principle  of  operation  of  RTD:  (a)  the  quantized  states  in  the  quantum  well  are  not  , 
aligned  with  occupied  states  in  the  emitter,  thus  no  current  flows,  (b)  the  RTD  is  “on  resonance”, 
electrons  can  go  from  the  emitter  to  the  collector  through  a  quantized  state 


itself,  utilizing  the  quantum  properties  of  nanostructures.  Two  ftindamental  noise 
mechanisms  that  must  be  considered  on  the  nanoscale  are  thermal  noise  and  shot  noise. 

Shot  noise  is  associated  with  the  fluctuations  of  emission  times  of  carriers  across 
a  limiting  barrier  in  any  electronic  device.  This  type  of  noise  is  the  dominant  source  of 
noise  in  bipolar  transistors,  p-n  diodes  and  hetero-barrier  diodes.  The  power  spectral 
density  of  such  shot  noise  is  given  by  S/  =  2eyl,  where  1  is  the  average  current,  e  is  the 
magnitude  of  the  electronic  charge,  and  y  the  shot  noise  factor  [1].  For  y  =  1 .0,  the 
phenomenon  is  statistically  Poissonian  and  called  full  shot  noise.  However  significant 
deviations  from  the  full  shot  noise  occur  in  double-barrier  resonant  tunneling  diodes 
(RTDs).  Shot  noise  reduction  as  low  as  y  =  0.35  [2]  at  77  K  has  been  measured, 
generating  interest  in  using  resonant-tunneling  structures  for  their  low-noise 
characteristics. 

This  research  features  the  first  shot  noise  calculation  for  multiple  (>2)  quantum 
well  resonant  tunneling  structures  that  includes  electrostatic-feedback  and  self-consistent 
coupling  between  Schrddinger  and  Poisson  equations.  It  will  support  the  claim  that  shot 
noise  reduction  effect  is  caused  by  the  electrostatic  repulsion  endured  by  injected  carriers 
due  to  the  charge  stored  in  the  quantum  well(s).  The  thesis  is  organized  as  follows.  The 
second  section  presents  the  operating  principles  and  fabrication  processes  of  RTDs.  The 
third  section  develops  in  detail  the  methodology  enabling  the  simulation  of  the  current- 
voltage  characteristics  of  resonant  tunneling  diodes  from  fundamental  principles. 
Fluctuations  around  the  average  values,  and  particularly  shot  noise,  are  studied  in  the 
fourth  section.  In  the  fifth  section  the  simulation  results  are  analyzed  and  compared  to 
experimental  results  for  an  InGaAs/AlAs  single- well  RTD.  Subsequently  the  shot  noise 
reduction  effect  is  studied  for  multiple-well  (up  to  4)  RTDs,  suggesting  that  properly 
designed  multiple-well  RTDs  will  display  even  lower  shot  noise  properties  than  single¬ 
well  devices.  In  the  final  section  the  use  of  multiple-well  RTDs  for  low-noise  application 
is  envisioned. 
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Fig.  3  Typical  I-V  curve  of  a  double-barrier  RTD 


2.  Noise  suppression  in  resonant  tunneiing  devices 

Resonant  tunneling  diodes  are  fabricated  from  layers  of  two  different  IIW 
semiconductor  alloys,  such  as  the  pair  Ino.53Gao.47As/AlAs.  A  double-barrier  resonant 
tunneling  diode  is  made  by  embedding  two  layers  of  the  material  with  the  larger  band  gap 
(AlAs)  within  layers  of  the  other  material  (InGaAs),  creating  a  quasibound  potential  well 
where  electrons  can  reside.  RTDs  are  made  with  wells  typically  4  to  6  nanometers  in 
width.  Whenever  electrons  are  confined  between  two  such  closely  spaced  barriers, 
quantum  mechanics  restricts  their  energies  to  one  of  a  discrete  set  of  quantized  levels. 

The  only  way  for  electrons  to  pass  through  the  device  is  to  "tunnel,"  quantum 
mechanically,  through  the  two  barriers.  The  probability  that  the  electrons  can  tunnel  is 
dependent  on  their  energy  relative  to  the  discrete  levels  in  the  well  of  the  device.  As 
illustrated  in  Fig.  2(a),  if  the  energy  of  the  incoming  electron  differs  from  the  levels 
allowed  inside  the  potential  well,  then  current  does  not  flow.  However,  when  the  energy 
of  the  incoming  electron  aligns  with  one  of  the  levels,  as  shown  in  Fig.  2(b),  the  electron 
is  said  to  be  "in  resonance"  with  the  well  and  current  flows  through  the  device. 

The  resulting  current-voltage  curve  (Fig.3)  features  a  negative  differential 
resistance  region  at  room  temperature,  making  the  RTD  a  practical  active  device.  P-N 
(Esaki)  tunnel  diodes,  which  display  the  same  kind  of  characteristics,  are  used 
commercially  in  microwave  oscillators,  mixers,  detectors,  and  multipliers.  RTDs  are 
displacing  Esaki  diodes  in  some  of  these  applications  because  of  their  superior  speed. 

2.1.  Fabrication  Procedure 

The  most  critical  part  of  the  fabrication  procedure  for  resonant  tunneling  diodes  is 
the  growth  of  the  wafer.  The  peak  current  density  of  RTDs  is  inverse-exponentially 
proportional  to  the  width  of  the  barriers.  So  the  thickness  of  the  grown  barriers  must  be 
carefully  monitored  (using  techniques  such  as  reflection  high  energy  electron  diffraction 
or  RHEED  that  provide  feedback  during  epitaxial  growth)  in  monolayer  accuracy.  Other 
requirements  are  to  avoid  the  diffusion  of  the  dopants  into  the  barriers,  to  smooth  the 
heterointerface  between  the  two  materials,  and  to  accommodate  the  residual  stress 
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Fig.  4  Left:  mesa  defining  the  RTD  —  Right:  zoom  on  the  under-etch  due  to  RIE 


between  the  two  semiconductor  alloys  (for  example  the  material  pair  Ino.53Gao.47 As/ A1  As 
exhibits  a  3.5%  lattice  mismatch). 

Once  the  RTD  wafer  is  grown,  fabrication  of  diodes  can  be  achieved  via  a  simple 
one-mask  process.  Grids  of  ohmic  metal  (e.g.  Ni/Ge/Au)  top  contacts  are  defined  by  a 
standard  photolithographic  liftoff.  The  metal  contacts  are  then  used  to  mask  the  sample 
during  the  etching  step  (e.g.  reactive  ion  etching  with  CCI2F2),  creating  mesa  devices 
(Fig.  4).  The  etching  process  is  stopped  about  100  nm  below  the  bottom  barrier  of  the 
RTD.  The  second  ohmic  contact  is  made  on  the  backside  of  the  wafer  or  around  the  mesa. 
If  probing  is  too  tricky  due  to  the  small  size  of  eontacts  or  if  the  devices  need  to  be 
packaged  for  RF  measurements,  additional  steps  are  required.  First,  a  dielectric 
deposition  must  be  made  by  PECVD  to  fill  the  volume  around  and  on  top  of  the  mesas, 
and  then  RIE  is  used  to  etch  openings  on  top  of  the  mesas.  Finally  metal  is  deposited  and 
patterned  to  define  large  landing  pads. 

2.2  Modeling  of  current-voltage  characteristics  of  RTDs 

Various  formalisms  of  quantum  transport  in  resonant  tunneling  diodes  have  been 
developed  over  the  years  [3].  In  this  work  the  goal  was  not  to  develop  an  exhaustive 
model  taking  every  physical  effect  into  account,  but  rather  to  develop  a  framework 
suitable  for  the  study  of  noise  properties  in  structures  with  thin  barriers.  In  sueh  structures 
the  electron  transport  is  largely  ballistic,  which  simplifies  the  analysis  significantly.  From 
a  first-principles  point  of  view,  the  RTD  is  an  open  quantum  system  in  which  electron 
transport  is  three  dimensional,  time  dependent,  time  irreversible,  dissipative,  and  many- 
body  interactive.  Both  particle  and  energy  exchanges  occur  frequently  with  the  outlying 
semiconduetor  material.  The  device  is  therefore  very  different  from  a  simple  isolated 
quantum  system,  where  a  conservative  Hamiltonian  and  the  boundary  conditions  for  the 
Schrodinger  equation  may  be  readily  formulated.  On  the  other  hand,  a  full  many-body 
formalism  may  not  be  necessary  for  calculations  of  specific  device  properties  since 
approximations  and  simplifications  can  be  used,  provided  that  the  major  transport 
processes  can  be  properly  modeled. 


100  nm  Ino.53Gao.47As  Nd=1.3  x  10^^  cm‘^ 

100  nm  Ino.53Gao.47As  Nd“2.0  x  10  cm 

5  nm  lno.53Gao.47As  Nd=2.0  x  lO'^cm*^ 

1 .4  nm  AlAs  UNDOPED 

6.1  nm  Ino.53Gao  47As  UNDOPED 

1.4  nm  AlAs  UNDOPED 

75  nm  Ino.53Gao.47As  Nd~1.0  x  10^  cm 

200  nm  Ino.53Gao.47As  Nd“2.0  x  10^^  cm  ^ 

200  nm  Ino.53Gao.47As  Np^l.S  x  10  cm 

InP  SUBSTRATE 

Fig.  5  Epitaxial-layer  profile  of  the  test  structure 

The  model  used  in  this  work  is  based  on  the  envelope  function  description  of 
electron  states.  Since  quantum  phenomena  must  be  included,  the  effective  mass 
Schrodinger  equation  is  invoked  to  calculate  the  electron  wave  function,  neglecting  the 
rapidly  varying  periodic  component.  The  primary  underlying  assumption  is  that  the 
transport  is  coherent,  meaning  that  there  is  conservation  of  the  transverse  crystal 
momentum  of  the  tunneling  carriers  and  that  phase-breaking  processes  are  neglected  in 
the  quantum-well(s)  region.  This  assumption  remains  reasonable  as  long  as  the  length  of 
the  tunneling  region,  and  especially  the  thickness  of  the  barriers,  is  negligible  in 
comparison  with  the  mean  free  path  of  electrons  under  thermal  excitation,  which  is  the 
case  for  the  structures  studied  in  this  work. 

All  through  this  section,  the  presentation  of  the  modeling  will  be  illustrated  by  a 
realistic  test  structure  whose  epitaxial-layer  profile  is  shown  in  Fig.  5.  A  few  remarks  can 
be  made  about  the  design  of  this  structure.  As  in  every  RTD,  both  ends  of  the  device  are 
heavily  doped  to  obtain  the  least  resistive  ohmic  contact  as  possible.  The  AlAs  barrier 
layers  and  the  Ino.53Gao.47As  quantum  well  are  undoped  and  surrounded  by  lowly  doped 
buffer  layers  to  preserve  the  purity  and  the  quality  of  the  barriers  and  heterointerfaces. 

The  buffer  layer  on  the  emitter  side  is  made  short  to  obtain  a  high  peak  current  density. 
The  collector  buffer  layer  is  made  long  because  the  capacitance  of  the  device  decreases 
with  the  thickness  of  this  layer  [6],  making  the  device  faster.  Oscillations  up  to  200  GHz 
have  been  obtained  with  this  structure  [4]. 

2.2.1  Material  parameters  and  energy-band  structure 

Accurate  material  and  energy-band  parameters  are  required  to  be  able  to  simulate 
quantitatively  any  nanostructure.  For  example,  because  the  peak  current  density  of  a  RTD 
is  exponentially  proportional  to  the  width  of  the  barriers,  a  good  knowledge  of  the  lattice 
constant  is  important  because  the  width  of  the  barriers  is  necessarily  a  multiple  of  the 
monolayer  length  of  the  material.  The  Ino.53Gao.47As/AlAs  material  system  presents  an 
additional  difficulty  because  these  materials  are  not  lattice-matched.  Because  the  AlAs 
barriers  are  very  thin,  it  is  assumed  that  all  the  stress  due  to  the  3.5%  lattice  mismatch  is 
endured  by  AlAs  [5].  Because  Ino.53Gao.47As  has  a  larger  lattice  constant,  the  lateral  stress 


7 


bond  struciure  AlAs 


Fig.  6  Band  structure  of  AlAs  -  The  dashed  curve  is  the  parabolic  band  structure,  the  plain 
curve  is  the  corrected  band  structure  used  in  the  simulation  —  The  left  part  is  the  complex  band 
structure,  the  right  part  is  the  real  one  -  The  zero  energy  is  the  edge  of  the  conduction  band 


is  tensional.  Thus  by  Hooke’s  law,  the  longitudinal  stress  is  compressive.  It  follows  that 
the  AlAs  band  gap  is  reduced  along  with  its  longitudinal  lattice  constant  and  conduction- 
band  offset.  The  parameter  values  used  in  the  simulation  are  listed  in  Appendix  1 . 

In  this  work,  only  electron  tunneling  through  the  T  valley  is  considered,  consistent 
with  the  fact  that  intersubband  scattering  is  neglected.  This  is  perhaps  the  most  critical 
simplification  in  comparison  with  more  rigorous  models,  but  it  is  also  what  makes  our 
model  elegant  and  easy  on  computer  resources.  This  is  often  called  the  single-band  model 
since  only  T-point  conduction-band  component  of  the  envelope  function  is  used  in  the 
tunneling  computations.  The  T-point  valence-band  components  are  assumed  to  be 
negligible. 

Using  a  single-band  model,  the  introduction  of  the  non-parabolicity  in  the  band 
structure  is  critical  to  obtaining  a  good  agreement  with  experiments  (Fig.  6).  The 
conduction  band  offset  for  Ino.53Gao.47As/AlAs  heterostructures  is  estimated  at  1 .2  eV.  So 
the  tunneling  of  electrons  through  the  AlAs  barriers  happens  very  deep  below  the  AlAs 
conduction  band  edge.  A  parabolic  imaginary  band  structure  will  have  the  effect  of 
strongly  underestimating  the  current  density  through  the  device  [7].  So  the  parabolic 
dispersion  relationship  is  generalized  to  an  imaginary  band  structure  that  connects  the 
conduction  band  and  the  light-hole  valence  band  [8]: 


E  =  E,- 


m 


r 

1  1 

V 

^  J 

with  rj  = 


1  1 
-T^ - 

m, 


m 


■Ih  J 


In  the  material  composing  the  quantum  well,  non-parabolicity  is  also  important  because  it 
has  a  great  impact  on  the  calculated  position  of  the  energy  levels  in  the  quantum  well. 
Non-parabolicity  is  introduced  via  an  energy  dependent  effective  mass  [9]: 


8 


Emitter  Region  Quantum  Well  Collector  Region 

Region 

Fig.  7  Typical  potential  profile  of  a  double-barrier  RTD  under  bias 


.  , .  .  /  .  f  m*(0) 

m  (E)  =  m  (0)(l  +  aE)  with  a=  1 - — 

V  ^0 

where  nio  is  the  vacuum  mass  of  an  electron  and  Eg  is  the  gap  energy.  These  corrections 
to  the  parabolic  dispersion  relationship  allow  one  to  avoid  the  use  of  a  multi-band  tight- 
binding  model,  and  appear  to  be  accurate  enough  to  simulate  the  main  physical  features 
of  resonant  tunneling  diodes. 

2.2.2  Computation  of  the  charge  and  current  density 

A  first  look  at  the  RTD  typical  potential  profile  (Fig.  7)  brings  out  the  general 
features  of  an  electron  wave  function  as  it  traverses  the  double-barrier  structure.  The 
highly  doped  emitter  acts  as  a  reservoir,  or  Fermi  “sea”,  of  electrons  for  tunneling 
through  the  potential  barriers.  The  electron  vmdergoes  multiple  reflections  between  the 
barriers  due  to  its  wave  nature,  leading  to  a  resonance  condition  of  the  tunneling  when 
approximately  one-half  of  a  de  Broglie  wavelength  equals  the  well  width.  The  high 
electric  bias  field  into  the  collector,  which  is  mainly  depleted  of  negative  charges,  then 
sweeps  the  electrons.  Also,  there  are  dimensional  changes  across  the  interfaces,  i.e.,  the 
three-dimensional  density-of-states  in  the  emitter  and  collector  versus  the  two- 
dimensional  density-of-states  in  the  quantum  well.  Thus,  in  modeling  the  resonant 
tunneling  diode,  three  regions  are  distinguished. 

The  emitter  and  the  collector  are  treated  the  same  way.  In  both  regions  it  is 
assumed  that  a  quasi-equilibrium  carrier  distribution  is  maintained  even  when  the 
structure  is  under  voltage  bias.  It  can  be  justified  by  the  fact  that  potential  barriers  tend  to 
isolate  these  regions  from  one  another  and  that  the  coupling  with  the  quantum  well  is 
weak,  even  on  resonance,  because  of  the  small  width  of  the  longitudinal  component  of 
the  density-of-states  in  the  quantum  well.  In  each  region,  the  quasi-Fermi  levels,  Ep  and 
Ef,  are  defined  by  the  requirement  for  space-charge  neutrality  over  the  entire  extent  of 
the  device.  The  Thomas-Fermi  local  equilibrium  treatment  of  the  conduction  electrons 
has  proven  robust  in  much  of  nanostructure  modeling  [10]  and  allows  the  calculation  of 


the  mobile  charge  in  the  emitter  and  the  collector; 


n(x)  =  iV,-^F, 

ylTT  ' 


k,T 


with 


yjx  dx 


+  e 


x-n 


N=2 


2nm  kgT 


However  this  approach  is  not  suitable  for  the  quantum  well  region  as  it  assumes  electrons 
are  distributed  as  a  3D  ftee  electron  gas. 

In  the  quantum  well  region,  electron  states  are  treated  purely  quantum 
mechanically.  The  one-dimensional,  time-independent,  single-particle,  single-band, 
effective  mass  Schrodinger  equation  is  used  in  the  direction  perpendicular  to  the 
quantum- well  layers: 


m  (x) 


V. 


+ 


V(x) 


where  E  is  the  longitudinal  energy  of  the  particle  considered,  Ec  is  the  bottom  energy  of 
the  conduction  band,  m* (z)  is  the  effective  mass,  and  V(z)  is  the  potential-energy  term 
including  the  effects  of  the  device  conduction-band  offset,  the  bias  voltage  applied  across 
the  device,  and  the  mobile  electronic  charge.  The  crucial  assumption  made  by  using  the 
one-dimensional  Schrodinger  equation  is  that  longitudinal  and  transverse  components  of 
the  energy  of  electrons  are  decoupled.  Both  the  longitudinal  energy  of  electrons  and  the 
length  (x-axis)  of  the  device  are  numerically  discretized.  Within  each  spatial  increment, 
the  mass  and  the  potential  are  assumed  to  be  constant  so  that  the  Schrbdinger  equation 
can  be  solved  with  piecewise-continuous  plane  waves  [1 1].  For  electrons  incident  from 
left  (emitter  region)  to  right  (collector  region),  we  have; 


^'~''[E,  x)  =  a'~''  (E,  x)  +  b'~'' (E,  x)  in  the  quantum  well  region, 

.  in  the  emitter  right  next  to  the  first  interface,  and 

in  the  collector  right  next  to  the  last  interface. 

As  indicated  above,  solutions  are  matched  to  plane-wave  functions  in  the  emitter  and  the 
collector  with  the  amplitude  of  the  incident  plane- wave  function  set  to  unity.  This 
solution  does  not  account  for  electron  flowing  from  the  collector  to  the  emitter.  For  such 
electrons: 
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x)  =  a'^~'  {£',  x)  (£',  x)  e''‘’^  in  the  quantum  well  region, 

« in  the  collector  right  next  to  the  last  interface,  and 
'Fg"^  =  f~'(^£')e~''‘‘’'‘'  in  the  emitter  right  next  to  the  first  interface. 

The  continuity  of  the  wave  function  'F^”'^(£,x)  and  of  its  derivative 

|^/w*(x  )]  d^'  x)/cbc  along  the  x-axis  gives  relationships  between  the  amplitude 
coefficients  of  the  wave  function: 


tVi 

where  the  subscript  j  refers  to  the  j  spatial  increment.  So  the  Schrodinger  equation  is 
solved  recursively,  from  right  to  left  (left  to  right)  for  electrons  flowing  from  emitter  to 
collector  (collector  to  emitter),  and  normalized  so  that  the  amplitude  of  the  incident  plane 
wave  is  unity.  The  statistical  properties  of  carriers  will  be  addressed  in  the  expressions 
used  to  calculate  the  charge  stored  in  the  well  and  the  current  density  across  the  structure. 
A  quantity  of  interest  for  the  current  density  calculation  is  the  transmission  probability  of 
an  electron,  i.e.  its  probability  to  cross  through  the  quantum-well  region  as  a  function  of 
its  initial  longitudinal  energy.  From  time  reversal  symmetry  and  the  probability-current 
continuity  principles  implied  by  the  quantum  mechanical  treatment,  we  obtain  useful 
relations  among  the  coefficients  C,  and  [12]: 
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r'-'(£)|  =  |r'-'(£)l 

4.r(£)i=*.r(^)i 

t|r'(£)f+lr'-'(£)f=i 

Thus  the  transmission  probability  is  the  same  in  both  directions  given  that  the  energy 
reference  is  the  same: 

T{E)  =  1  -p-''  (£)|^  =  1  -|r'"'  (£)|^ 

An  important  remark  about  the  discretization  of  the  Schrodinger  equation  has  to  be 
made.  Although  the  spatial  increments  are  made  uniform  for  technical  reasons  explained 
later,  the  energy  discretization  is  made  “adaptive”.  The  transmission  function  as  well  as 
the  wave  function  feature  sharp  peaks  around  the  values  of  the  energy  levels  in  the 
quantum  well(s).  The  energy  increments  are  made  finer  around  those  energies  to  achieve 
a  good  resolution  of  the  transmission  peak.  This  is  done  by  linking  the  energy  increments 
to  the  slope  of  the  transmission  function.  This  adaptive  stepping  technique  saves  a 
tremendous  amount  of  computation  time  and  refines  the  calculation  aroimd  the  energies 
that  contribute  most  to  the  current  density  and  the  electronic  charge  density. 

The  current  density  can  be  evaluated  using  the  Tsu-Esaki  formula  [13],  which  can  be 
explained  as  follow  [14].  Planes  waves  incident  from  the  emitter  having  longitudinal 
energy  E  and  transverse  energy  Et  carry  a  statistical  'weight  fe(E+E,),  where _/e  is  the 
Fermi-Dirac  distribution  in  the  emitter.  They  give  rise  to  an  incident  electron  current 
contribution  from  a  unit  volume  of  k-space: 


dE 


with  the  factor  2  in  the  numerator  accounting  for  spin  and  the  (27i)^  for  the  density  of 
states  in  k-space.  The  last  term  represents  the  group  velocity  of  the  plane  wave.  To 
contribute  to  the  transmitted  flux  this  plane  wave  must  not  only  couple  to  a  state  on  the 
other  side  of  the  structure,  but  because  electrons  are  fermions,  that  state  must  be  empty. 
The  transmitted  electron  current  contribution  is  then: 


2e 


-/.(£  +  £, )[l-/,(£  +  AK  +  £,)]r(£)i 


^  dE^ 


\^Kj 


where  fc(E+AV+E,)  is  the  Fermi-Dirac  distribution  in  the  collector  and  A  V  the  potential 
energy  drop  between  the  emitter  and  the  collector  or  the  longitudinal  kinetic  energy 
gained  by  an  electron  after  crossing  the  quantum- well  region.  If  all  energy  levels  are 
referenced  to  the  same  point,  we  have: 
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El-Ei=V, 


bias 


f,(E*hV  +  E,)  = 


1  +  exp 


1  +  exp 


e+e,-e; 

k,T 


e+e,-e; 


n-l 


k,T 


Vbias  is  the  fraction  of  applied  voltage  that  drops  across  the  structure.  The  total  electron 
current  density  flowing  from  left  to  right  is  obtained  by  integrating  over  the  longitudinal 

and  transverse  wave  vector,  kg  and  k, ,  of  electrons  in  the  emitter. 


dk^dk, 


To  be  able  to  partially  integrate  this  equation  analytically,  we  assume  that  the 
transverse  dispersion  relationship  is  parabolic: 


2 1,  2 


E.  = 


rk, 

2ml 


After  a  change  of  variable,  the  expression  used  in  the  simulation  is  obtained: 


W-r 


J7t£)ln 

(  E^-E 

1  +  e 

E]r-E 
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dk^ 


l-e 
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The  same  kind  of  analysis  allow  one  to  obtain  the  current  density  for  electrons  flowing 
from  right  to  left: 

BE' 

{  ^F-tL  \  1 

j,.,  ^em%T 


271  E 


2*3 


/Jt£)ln 

l  +  e 

Ef-E 

[l  +  e'‘‘^  J 

dr 


dr 


l-e 


The  total  electron  current  density  is  given  by: 

It  can  be  checked  that  this  current  density  goes  to  zero  when  no  bias  voltage  is  applied. 
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and  that  the  component  becomes  negligible  as  an  increasing  positive  bias  voltage  is 
applied. 

The  electron  density  in  the  quantum  well  region  is  also  calculated  by  considering  two 
streams  of  electrons  impinging  jfrom  the  emitter  and  the  collector,  respectively  [15].  For 
the  left-to-right  electron  flow,  the  modulus  square  of  the  wave  function  is  used  to 
determine  how  electrons  in  the  Quasi-Fermi  sea  on  the  emitter  side  spread  into  the 
quantum  well  region: 


(^)  =  +  dk^dk, 

(2;r)  ' 


J|'F'-^|'(£,x)ln 


\  +  e 


E’f-e\ 

k„T 


dK 


Again  the  transverse  component  of  the  wave  vector  has  been  integrated  analytically.  A 
similar  expression  is  derived  for  the  stream  of  electrons  flowing  from  right  to  left: 


|'F^-'p(£:,jc)ln 


(  E^r-E\ 


dK 


and  (x)  =  n'  "  (x)  +  n''  ‘  (x) 


It  is  important  to  remark  about  the  boundaries  of  integration  in  the  expressions  of 
the  current  density  and  the  electron  density.  When  the  RTD  is  biased  positively,  more  and 
more  charge  accumulates  just  before  the  first  heterointerface  in  the  low-doped  part  of  the 
emitter.  The  potential  profile  in  this  section  of  the  structure  presents  a  triangular-like 
quantum  well  feature,  and  quantized  levels  should  be  formed.  The  inclusion  of  this 
additional  quanti2:ation  is  tedious  and  not  compatible  with  the  Tsu-Esaki  formula  (which 
assumes  a  3D  Fermi  distribution  in  the  emitter)  [16]-[17].  Therefore  for  the  sake  of 
simplicity,  only  states  above  the  top  of  the  conduction-band-edge  maximum  are 
accounted  for  in  the  current  and  electron  density  calculations  (as  illustrated  in  Fig.  7). 
This  simplification  is  not  as  bad  as  it  first  appears.  Electrons  with  a  small  longitudinal 
energy  have  only  a  small  contribution  to  the  current  density,  as  their  group  velocity  is 
also  small.  We  believe  this  approximation  does  not  affect  much  the  overall  current 
density,  though  it  has  an  impact  on  the  voltage  range  of  the  negative  differential 
resistance  region. 

We  now  have  physical  models  for  the  current  density  and  the  electron  density 
along  the  entire  structure.  Both  have  to  be  coupled  to  the  Poisson  equation  to  determine 
the  physical  characteristics  of  the  resonant  tunneling  diodes. 

2.2.3  Computational  method 


The  model  presented  above  was  used  to  determine  the  electronic  charge  density  from 
the  potential  profile.  With  the  Poisson  equation  these  equations  form  a  self-consistent  set 
of  equations: 
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where  AEf  is  a  constant  related  to  the  difference  of  doping  between  the  ends  of  the 
collector  and  emitter. 

The  first  step  is  to  discretize  the  Poisson  equation  to  transform  the  differential 
equation  into  a  linear  equation  [18]: 
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«F(0)  =  0 


Ax 

F,=0 

V,=V^-K^, 


bias 


-  =  fj  0<j<N-\ 

+  AEp 


Because  of  the  uniform  spatial  increments,  the  Poisson  equation  can  then  be  expressed  in 
a  matrix  form: 


Ay  =  Ax^f-b 


'-2  1 
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•.  ••  1 
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v  = 

/v-J 

J  = 

>b  = 

The  matrix  A  can  be  easily  inverted  and  we  have  a  direct  linear  relationship  to  calculate 
the  potential  profile  from  the  charge  profile  along  the  structure  at  a  given  voltage  bias. 

We  now  have  a  system  of  two  algebraic  equations,  which  explicitly  relate  the 
charge  density  and  potential  profile  to  one  another,  to  resolve.  The  simplest  convergence 
algorithm  one  can  think  of  is  presented  in  Fig.  8  [19]-[20].  For  a  given  bias  voltage,  a 
first  guess  of  the  potential  profile  is  made.  The  charge  density  along  the  structure  is 
calculated.  Then  a  new  potential  profile  is  calculated  thanks  to  the  Poisson  equation.  If  all 
the  points  of  this  new  potential  profile  are  within  a  pre-determined  margin  of  the  old 
potential  profile,  the  system  has  converged.  If  it  is  not  the  case,  this  new  potential  is 
launched  back  into  the  loop  until  convergence  is  achieved.  Once  a  solution  is  obtained  for 
a  given  voltage  bias,  this  solution  makes  a  good  first  guess  for  the  following  bias  point. 

Actually  this  algorithm  does  not  give  satisfying  results.  It  is  too  unstable  and 
solutions  tend  to  diverge.  Even  the  addition  of  a  relaxation  factor  is  not  sufficient  to 
stabilize  the  system.  That  is  the  reason  why  the  addition  of  an  acceleration  method  [21]- 
[22]  has  proven  to  be  necessary.  The  whole  process  of  solving  the  Poisson  and  charge 
density  equations  may  be  regarded  as  an  operator  //transforming  F*  into  W^: 


If*  =  //f  * 


where  f*  is  the  k‘'’  iteration  of  the  unknown  potential  profile  and  If*  is  the  potential 
profile  calculated  from  f*.  The  basic  algorithm,  where  no  acceleration  method  is  present, 
can  be  expressed  as  follow: 

_  jyk 

The  idea  underlying  the  acceleration  method  is  to  use  If*,  along  with  the  V’s  and  W’s 
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stored  during  the  previous  iteration,  to  evaluate  the  new,  optimized,  iteration  as: 

The  purpose  of  the  function  F  is  to  build  up  an  iterative  sequence  that  converges  more 
rapidly  than  the  basic  sequence.  The  acceleration  method  works  as  follows.  Define 


7=1 

where  Mis  the  number  of  previously  stored  iterates,  and  a/’s  are  free  parameters.  The 
latter  are  chosen  in  order  to  minimize  the  residual 


^  /=! 


2 


where  A'^is  the  dimension  of  f*.  In  turn,  minimizing  requires  the  algebraic  system  of 
order  M 

dR” 


0a; 


=  0  y=i,...,M 


to  be  solved.  Finally  the  new  iterate  is  defined  as 


=«*  +  /?  (6* -a*)  O<y0<l 

where  p  is  the  relaxation  parameter.  The  simulations  have  been  carried  out  with  M=  5 
and  P  has  been  varied  according  to  the  fastest  convergence  time  for  each  structure  under 
investigation. 

The  convergence  criterion  used  is  that  the  calculated  potential  must  be  within  0.5 
meV  of  the  previous  calculated  potential.  Self-consistent  solutions  are  obtained  in  a  few 
tens  of  second  on  a  2  GHz  machine  after  5  to  20  iterations,  depending  on  the  length  of  the 
structure  and  the  bias  voltage  at  which  the  calculation  is  made.  The  modeling  of  the 
current-voltage  characteristics  of  resonant  tunneling  diodes  is  only  the  first  step  of  the 
simulation  procedure  for  determining  the  noise  characteristics,  presented  in  the  next 
section. 

2.3  Noise  mechanisms  in  RTDs 

Two  noise  mechanisms  must  be  considered  in  resonant  tunneling  diodes:  thermal  noise 
and  shot  noise. 
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Fig.  9  Band  bending  model  of  single-well  structure 


2.3.1  Thermal  noise 

Thermal  noise  is  caused  by  the  scattering  of  carriers  with  a  lattice  and  is  directly 
proportional  to  the  temperature  of  the  thermodynamic  bath  in  which  they  reside.  Once  the 
I-V  characteristic  of  the  RTD  is  known,  the  spectral  density  of  current  fluctuations  for 
thermal  noise  can  be  immediately  estimated  from  Nyquit’s  generalized  formula; 

Sj  =  4ksTG 

where  ke  is  Boltzman’s  constant,  T  is  the  device  temperature,  G  -  Re{Y}  and  Y is  the 
differential  admittance.  Thus  it  appears  that  the  thermal  noise  power  increases  with  the 
absolute  value  of  the  slope  of  the  I-V  curve. 


2.3.2  Shot  noise 

Shot  noise  is  a  type  of  fluctuation  in  the  electrical  current  that  is  observed  in  many 
electronic  devices  operating  out  of  thermodynamic  equilibrium.  In  most  cases  it  is 
associated  with  randomness  in  the  flux  of  carriers  crossing  a  barrier  in  the  device.  The 
source  of  this  randomness  is  thermal  fluctuations  in  the  rate  of  carrier  injection  into  the 
barrier.  At  frequencies  small  compared  to  the  reciprocal  of  the  transit  time  across  the 
active  layer,  the  power  spectral  density  of  such  shot  noise  is  given  by  Si  -  2eyl,  where  1  is 
the  average  current,  e  is  the  magnitude  of  the  electronic  charge,  and  y  is  the  shot  noise 
factor. 

In  the  majority  of  devices  it  is  found  that  y  =  1.0,  which  is  the  condition  of 
normal,  or  full,  shot  noise.  Some  examples  are  p-n  diodes,  Schottky  diodes,  and  single 
hetero-barrier  diodes.  In  these  devices  the  potential  energy  barrier  has  transmission 
characteristics  that  are  independent  of  the  current  magnitude.  The  carrier-injection  times 
into  the  active  layer  obey  a  Poissonian  distribution.  This  implies  that  these  injection  times 
are  statistically  independent  and  that  the  transit  of  a  carrier  across  the  active  layer  is  an 
event  fully  uncorrelated  to  the  transit  of  other  carriers.  Therefore,  the  arrival  times  of 
carriers  on  the  collector  side  of  the  active  layer  obey  Poisson  statistics,  leading  to  the  full 
shot  noise  in  the  terminal  current  fluctuations. 

In  devices  having  a  potential  barrier  whose  transmission  depends  on  the  current, 
significant  deviation  from  frill  shot  noise  can  occur.  The  classic  example  is  the  space- 
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charge-limited  vacuum  diode  in  which  the  barrier  consists  of  a  spatial  distribution  of 
electrons  located  near  the  cathode  [23].  In  double-barrier  resonant  tunneling  diodes,  shot 
noise  reduction  (y  <  1.0)  [24]-[27]  and  shot  noise  enhancement  (y  >  1.0)  [27]  have  been 
observed.  The  primary  assumption  of  the  following  derivation  is  that  the  physical  basis  of 
the  shot  noise  reduction  and  enhancement  is  the  modulation  of  the  double-barrier 
transmission  probability  by  the  charge  stored  in  the  quantum  well.  Notably  the  charge  in 
the  well  induces  correlations  between  the  fluctuations  of  the  incident  current  in  different 
energy  intervals.  This  effect  will  be  more  clearly  explained  in  the  presentation  of  the 
analytical  model  of  shot  noise  used  in  this  work. 

2.4.  Modeling  of  shot  noise  in  RTDs 

2.4.1.  Analytic  Modeling 

Let  us  consider  a  single- well  resonant  tunneling  diode  biased  with  a  large  enough 
voltage  so  that  a  wide  depleted  space-charge  region  forms  on  the  collector  side  and 
excess  electrons  accumulate  on  the  emitter  side.  Under  this  bias  condition,  the  current 
through  the  diode  consists  of  electrons  that  are  incident  on  the  double-barrier  structure 
from  the  emitter  side  and  are  transmitted  through  the  structure  by  resonant  tunneling 
through  the  first  or  higher  quasi-bound  levels  of  the  quantum  well.  Ji  denotes  the  total 
time-averaged  electronic  flux  incident  on  the  double-barrier  structure  from  the  emitter 
side,  and  Jt  denotes  the  transmitted  flux.  These  fluxes  can  be  decomposed  into 
components  J„'  and  J„' contained  within  specified  intervals  of  energy  /IE„  along  the 
tunneling  direction  (Fig.  9).  The  energy  is  referenced  to  the  conduction  band  edge  on  the 
emitter  side.  The  magnitude  of  the  time-averaged  electrical  current  through  the  diode  is 
related  to  the  incident  flux  by: 

n  n,nj  m  ni 


where  is  the  transmission  probability  through  the  double-barrier  structure  of  an 
electron  incident  in  the  mth.  energy  interval  and  transmitted  to  the  nth  interval,  and  Ij 
is  the  electrical  current  flowing  in  the  mth  interval  on  the  emitter  side.  We  assume  the 
resonant  tunneling  is  elastic,  so  that 

r„„,  =  0  for  m  12  and  J'„  =  T„J'„ 

However  the  formalism  developed  is  valid  independent  of  the  coherent  or 
sequential  nature  of  the  resonant  tunneling.  For  the  purpose  of  analyzing  the  shot  noise, 
fluctuations  in  the  current  and  the  fluxes  are  defined  in  the  usual  sense  as  the  difference 
between  the  instantaneous  and  the  average  value  of  each  quantity.  All  fluctuations  are 
assumed  to  occur  under  the  condition  of  an  ac  short  circuit  across  the  diode. 

The  source  of  the  shot  noise  is  the  fluctuation  of  the  incident  flux  in  each  energy 
interval.  A  given  component /1J„'  has  two  effects  on  the  transmitted  current.  First,  it 
contributes  directly  to  the  transmitted  current  through  the  term  Second  it 

modulates  the  transmission  probability  in  other  energy  intervals  by  varying  the  charge 
density,  and  hence  the  potential  profile,  in  the  double-barrier  structure.  To  first  order,  the 
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fluctuation  in  the  diode  current  due  to  fluctuations  in  all  intervals  is 


n  n,m  m 


Since  the  fluctuation  AJ„'  occurs  independently  of  T„,  one  can  define  A/®  =  .  It 

follows  that 

n 

where  =  J‘„  {dTjdJl^T^,^ .  The  quantity  C„m  is  an  element  of  a  matrix  that  represents 

the  modulation  of  the  transmission  at  each  energy  by  fluctuations  of  the  incident  current 
at  other  energies. 

The  low-frequency  spectral  density  of  the  total  current  fluctuation  is  found  to  b 


W-_L 

Sf  Sf 


V  » 


2  A/;  +xa/;Ic„a/:+ 


This  expression  is  simplified  with  the  following  assumptions.  First,  incident  current 
fluctuations  in  different  longitudinal  energy  intervals  are  imcorrelated  [23] 


ai‘ai‘ =s  Fa/'T 

n  m  nm  L  «  J 

where  d„m  is  the  Kronecker  delta  function.  Second,  the  incident  current  component 
fluctuations  are  Poisson-distributed  and  are  thus  related  to  the  average  component  by  the 
well-known  Schottky  relation  [1]: 


[a/.*]  =2ei:sf 

The  application  of  these  relations  leads  to  the  result 


Sj  =  lei 


i+2r'Xc„/:+E  Ic, 


nm 
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In  the  limit  that  AE„  — >  0  for  all  n,  y  can  be  written  in  integral  form  as 


Y 


-MI-'  l[\/C(E,B)]dI\EYl-'  £,[ \dC(E,E')]dI\E') 


Where  dC{E.,E')=  dJ  \E)[dT  {E)/ dJ \E')\t ''  (£') 
and  1=  [  dr  (E)  =  [  eT(E)  dJ\E) 
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Regular  Diode  Resonant  Tunneling  Diode 


Fig.  10  In  an  RTD  the  charge  stored  in  the  quantum  well  is  responsible  for  the  correlations 
between  current  fluctuations  at  different  energies,  leading  to  deviation  from  the  full  shot  noise. 


This  equation  is  the  basis  of  the  following  simulations  to  investigate  shot  noise  in 
multiple- well  resonant  tunneling  diodes. 

2.4.2.  Numerical  Modeling 

The  numerical  calculation  of  the  shot  noise  factor  requires  the  calculation  of  each 
element  of  the  correlation  matrix  dC{E,  E') .  The  only  part  that  is  not  readily  available  is 

the  derivative  5r(£)/6l7'(£') .  This  term  represents  the  impact  of  an  incident  current 

fluctuation  at  energy  on  the  transmission  probability  at  energy  E,  which  is  related  to 
the  amount  of  charge  injected  into  the  well(s)  due  to  this  current  fluctuations  (Fig.  10).  It 
is  calculated  using  the  following  procedure  using  the  chain  rule  of  differential  calculus: 

dT{E)  ET{E)  _ET{E)  An{x) 
dJ‘{E')  ~  Aj{e')  ~  An(x)  AJ‘{E') 
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First  a  current  fluctuation  /SJ{E')  is  properly  picked  so  that  the  resulting  shot  noise  factor 
is  independent  of  this  choice.  This  is  accomplished  by  choosing  LJiE’)  so  that  the 
amount  of  charge  injected  due  to  this  current  fluctuation  is  kept  a  constant  fraction  of  the 
average  amount  charge  in  the  well  at  the  bias  voltage  at  which  the  computation  is  made. 
This  fluctuation  gives  rise  to  a  fluctuation  An(x)  of  the  charge  in  the  quantum  well.  By 
identification  with  the  current  density  formula,  we  have: 
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Fig.  12  Epitaxial-layer  profile  of  the  test  structure  made  from  Ino.53Gao.47As/AlAs.  Note  that  the 
thicknesses  of  the  top  barrier  and  quantum  well  are  slightly  less  than  the  nominal  (MBE  growth) 
values,  consistent  with  the  theoretical  results  given  here. 


and  A?7(x)  = 


dJ‘{E')  ^  ’ 


Here  we  assumed  that  only  current  flowing  from  emitter  to  collector  gives  rise  to  shot 
noise,  which  should  be  a  good  approximation  as  long  as  the  bias  voltage  is  not  too  close 


to  zero. 

The  charge  fluctuation  A«(x)  modifies  the  potential  profile  of  the  device.  If  we 
assume  the  depletion  region  in  the  collector  is  not  affected  by  this  charge  fluctuation 
(because  positive  charges  are  not  affected  by  the  fluctuation),  the  Poisson  equation  can  be 
integrated  using  as  boundary  conditions  the  value  of  the  electric  field  at  the  collector  / 
quantum  well  region  interface  and  the  value  of  the  potential  at  the  emitter  end  as 
boundary  conditions.  A  new  transmission  function  is  then  obtained  from  this  modified 
potential  profile,  allowing  the  calculation  of  Ar(£).  This  numerical  procedure  can  be  very 
intensive  in  computer  resources  because:  (1)  the  correlation  matrix  is  not  symmetric  and 
each  matrix  element  must  be  calculated  individually,  and  (2)  the  number  of  matrix 
elements  is  equal  to  the  square  of  the  number  of  energy  intervals  and  the  number  of 
energy  intervals  must  be  a  three-to-four  digit  number  to  obtain  a  reasonable  precision. 


2.4.3.  Physical  Interpretation 


A  typical  shot  noise  characteristic  of  a  single- well  RTD  is  presented  in  Fig.  11. 
These  data  [2]  were  measured  for  a  GaAs/AlGaAs  structure  at  77  K  and  are 
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region,  shot  noise  reduction  (y<1)  is  observed.  In  the  negative  differential  resistance 
region,  shot  noise  enhancement  (y>l)  is  observed.  Based  on  the  formalism  developed  it 
can  be  explained  as  follows.  A  positive  fluctuation  of  the  incident  current  injects  more 
mobile  electronic  charge  in  the  quantum  well.  Through  the  Poisson  equation,  this  injected 
charge  modifies  the  potential  profile  of  the  structure,  causing  the  quasi-bound  levels  of 
the  quantum  well  to  shift.  This  then  impacts  the  transmission  probability  of  electrons  in 
the  emitter  in  a  manner  depending  on  the  bias  region  of  the  device. 

In  the  positive  differential  resistance  region  (PDR),  a  positive  current  fluctuation 
shifts  up  the  quasi-bound  level  of  the  quantum  well  to  an  energy  at  which  there  are  fewer 
electrons  to  tunnel  from  the  emitter  side.  This  is  because  in  the  PDR  region  the  quasi¬ 
bound  level  is  aligned  with  the  Fermi  sea  on  the  emitter  side  so  an  increase  in  energy 
corresponds  to  a  reduction  in  the  Fermi-Dirac  probability.  Thus  the  instantaneous  total 
current  is  decreased  by  the  fluctuation,  leading  to  shot  noise  reduction. 

In  the  negative  differential  resistance  region  (NDR),  it  shifts  up  the  quasi-bound 
level  to  an  energy  at  which  there  are  more  electrons  to  tunnel.  This  is  because  in  the  NDR 
region  the  quasi-bound  level  is  aligned  below  the  conduction-band  edge  of  the  Fermi  sea 
on  the  emitter  side  so  an  increase  in  energy  corresponds  to  an  increase  in  Fermi-Dirac 
probability.  An  increase  of  the  instantaneous  total  current  is  then  caused  by  the 
fluctuation,  leading  shot  noise  enhancement. 

It  is  also  appropriate  to  ask  what  the  meaning  is  of  shot  noise  reduction  and 
enhancement  in  a  more  general  sense.  When  there  is  full  shot  noise  (y=1),  the  emission 
and  collection  times  of  carriers  through  the  structure  follow  a  Poissonian  distribution. 
Thus  they  are  both  fully  independent  statistical  variables.  It  means  the  collection  of 
carriers  is  a  fully  random  process  and  each  collection  event  is  independent  from  the 
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Fig.  14  Design  of  the  three  InGaAs/AlAs  RTDs  under  investigation 


others.  The  introduction  of  the  quantum  well  causes  correlations  between  the  collection 
times  of  carriers.  In  the  case  of  shot  noise  reduction,  these  correlations  are  negative  and 
tend  to  reduce  the  randomness  of  the  collection  of  carriers.  In  the  limit  that  the  shot  noise 
is  totally  suppressed,  the  carriers  would  be  collected  at  a  perfectly  constant  rate,  a  carrier 
being  collected  exactly  every  x  seconds  where  r  =  e/I  and  /  is  the  average  current 
through  the  structure.  Shot  noise  enhancement  is  more  difficult  to  interpret.  Since  normal 
shot  noise  is  already  a  fully  random  process,  shot  noise  enhancement  cannot  be  regarded 
as  a  more  random  process.  Let  us  just  say  that  the  correlations  between  the  collection 
times  of  carriers  become  positive,  and  thus  the  fluctuations  of  the  device  terminal  current 
are  amplified.  This  results  in  an  increase  of  the  noise  spectral  density  from  the  device. 


2.4.4.  Comparision  of  numerical  model  and  experiment 


The  previous  sections  developed  an  elaborate  model  to  study  the  noise  properties 
of  resonant  tunneling  diodes.  The  validity  of  the  model  must  be  confirmed  by  comparison 
of  simulation  to  experiments.  The  present  analysis  is  limited  to  coherent  tunneling  at 
room  temperature  because  of  the  relative  simplicity  of  the  coherent  case  and  because 
coherent  tunneling  generally  defines  the  best  possible  electrical  performance  of  RTDs. 
The  test  structure  used  to  validate  the  model  was  originally  used  to  fabricate  an  oscillator 
[4]  and  the  data  on  shot  noise  are  unpublished.  The  test  structure  is  an 
Ino.53Gao,47As/AlAs  double-barrier  RTD  at  T  =  300  K.  The  structure  was  grown  by 
molecular-beam  epitaxy  and  fabricated  into  4-pm-diam  mesa-isolated  diodes.  The 
material  structure  is  shown  in  Fig.  12.  The  barriers  were  made  very  thin  and  the  doping 
concentration  on  the  emitter  (top)  side  very  large  to  obtain  high  current  density,  as 
evidenced  by  the  peak  value  of  approximately  2.5x10^  A-cm'^  in  the  experimental  I-V 
curve  of  Fig.  13(a).  For  the  present  purposes,  the  noise  power  from  the  packaged  diodes 
was  measured  around  1  GHz.  The  thermal  noise  was  subtracted  from  the  total  noise 
power  to  obtain  the  shot  noise  factor  at  selected  bias  points  as  plotted  in  Fig.  13(b).  As  in 
previous  measurements,  the  shot  noise  is  suppressed  (y  <  1 )  at  bias  voltages  below  the 
NDR  region  and  is  enhanced  above  the  NDR  region.  In  the  present  measurements  no 
results  were  obtained  within  the  NDR  region  because  of  the  difficulty  in  suppressing 
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Fig.  1 5.  Illustration  of  the  design  guideline  for  multiple-well  RTD 


oscillations. 

The  predictions  of  the  present  simulations  are  plotted  against  experiment  in  Fig. 
13.  The  agreement  is  good  for  both  the  I-V  characteristic  and  shot  noise  factor.  The 
voltage  at  which  the  NDR  region  occurs  is  well  predicted  given  the  facts  that:  (1)  the 
plateau  observed  experimentally  between  1.9  V  and  2.5  V  is  probably  caused  by  self¬ 
rectification  of  oscillations  with  the  external  circuit  which  distort  the  NDR  region,  and  (2) 
as  pointed  out  in  reference  [16]  an  overestimation  of  the  current  peak  and  an  abrupt  NDR 
region  is  to  be  expected  when  the  assumption  of  the  decoupling  of  the  longitudinal  and 
transverse  energies  is  made.  Moreover  the  simulation  exhibits  some  bistability  effect  in 
the  NDR  region,  not  shown  in  Fig.  13  but  addressed  in  Appendix  2,  consistent  with  the 
experimental  results.  The  shot  noise  factor  is  also  qualitatively  and  quantitatively  well 
predicted.  Notably  the  minimum  measured  shot  noise  factor  is  y"”"  =  0.45  and  the 
minimum  simulated  value  is  y"”"  =  0.39. 

The  bulk  of  the  disparities  between  simulation  and  experiment  can  be  attributed 
to  different  causes.  First  the  physical  model  treats  the  emitter  and  the  collector  semi- 
classically,  and  the  quantum  well  region  fully  quantum  mechanically.  A  quantum 
mechanical  treatment  of  the  entire  structure  would  be  more  accurate  but  requires  the 
introduction  of  heuristic  scattering  potentials  to  obtain  charge  conservation,  and  raises 
issues  of  boundary  conditions  of  open  systems  that  are  still  not  fully  understood  today. 
The  main  shortcoming  of  the  present  treatment  is  that  the  quantum-size  effects  related  to 
the  formation  of  the  accumulation  layer  at  the  emitter  /  double-barrier  interface  are 
neglected.  We  believe  it  is  a  small  effect,  although  it  is  likely  to  change  the  position  of 
the  NDR  region.  Scattering  is  also  neglected  but  that  seems  to  be  a  reasonable 
assumption  for  tbe  structures  under  investigation  that  have  very  thin  potential  barriers. 
Finally  the  single-band  model  does  not  account  for  tunneling  through  the  X-valley.  But 
that,  too,  may  be  negligible  because  of  the  thin  barriers  [28]-[29]. 

Despite  its  many  simplifications,  our  model  displays  reasonable  agreement  with 
experiments,  and  thus  allows  us  to  study  shot  noise  in  more  complicated  (multiple-well) 
structures. 
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2.5.  Simulation  of  shot  noise  in  multiple-well  RTDs 

2.5.1.  InGaAs/AlAs  Structures 

The  shot  noise  factor  and  I-V  characteristics  are  computed  at  room  temperature 
for  three  different  Ino.53Gao.47As/AlAs  structures  all  sharing  the  same  emitter  and 
collector  layers  (Fig.  14)  and  having  potential  barriers  1.4  nm  wide.  Structure  A1  is  a 
double-barrier  single- well  RTD  with  a  6.1  nm  wide  well,  structure  A2  is  a  triple-barrier 
RTD  with  6.1  nm  and  4.5  nm  wells,  and  structure  A3  is  a  quadruple-barrier  RTD  with  6.1 
nm,  4.5  nm  and  3.6  nm  wells.  Structures  A2  and  A3  are  designed  with  graded  well  widths 
(the  narrowest  well  on  the  collector  end)  so  that  the  ground-state  energies  of  the  three 
wells  gradually  come  into  co-alignment  as  the  bias  voltage  is  increased.  With  this  so- 
called  graded  quantum- well  technique  [30]-[31]  (Fig.  15),  the  peak  current  density  can  be 
made  usefully  high  if  the  co-alignment  occurs  while  the  states  are  dropping  through  the 
Fermi  sea  on  the  emitter  side.  The  wells  and  barriers  of  all  three  structures  are  assumed  to 
be  very  thin  to  maximize  the  likelihood  of  coherent  transport  without  scattering  and, 
thereby,  justify  a  key  assumption  of  our  model. 

The  predictions  of  our  model  for  structures  Al,  A2,  and  A3  are  plotted  in  Figs. 

16, 17,  and  18,  respectively.  Each  structure  exhibits  about  the  same  peak  current  density 
of  around  2x10"^  A-cm’^.  The  shot  noise  factor  curves  are  qualitatively  similar,  but  the 
shot  noise  reduction  effect  is  greatly  enhanced  in  the  multiple- well  structures.  The 
minimum  shot  noise  factor  computed  in  the  PDR  region  of  (the  single-well)  structure  Al 
is  0.81.  This  y^'"  drops  to  0.44  and  0.24  in  structures  A2  and  A3,  respectively.  We 
attribute  this  decrease  to  the  fact  that,  in  structures  having  the  same  peak  current  density, 
the  total  amount  of  charge  in  the  quantum  well  region  increases  with  the  number  of  wells. 
This  increases  the  electrostatic  repulsion  encountered  by  the  injected  carriers  and  thus 
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increases  the  negative  correlations  responsible  for  the  shot  noise  reduction  effect. 


Fig.  19  shows  the  minimum  shot  noise  factor  versus  number  of  wells  in  the 
structure.  There  is  clearly  a  decreasing  trend  that  leads  us  to  ask  whether  this  result  could 
significantly  impact  the  overall  noise  performance  of  a  real  device  made  from  such 
structures.  Because  most  devices  do  not  operate  under  thermodynamic  equilibrium,  there 
is  no  fundamental  limit  (equivalent  to  Shannon  theorem)  for  the  shot  noise  reduction 
effect.  But  this  argument  neglects  the  contribution  from  thermal  noise. 
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Bias  Voltage  (V) 

Fig.  18  Structure  A3:  Triple- well  RTD  -  (a)  current  density  and  (b)  shot  noise  factor  vs.  bias 
voltage.  A  minimum  shot  noise  factor  of  0.24  is  achieved.  The  shot  noise  factor  curve  exhibits 
three  peaks  as  expected  at  0.67  V,  0.73  V  and  0.75  V,  The  small  undulations  in  the  shot  noise 
factor  are  caused  by  numerical  truncation. 
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Fig.  20  Shot  noise  power  spectral  density  (solid  curve)  and  corresponding  thermal-noise 
spectral  density  (dotted  curve)  plotted  for  (a)  structure  A1  (single- well  RTD)  and  (b) 
structure  A3  (triple-well  RTD). 


From  Figs.  16, 17,  and  18,  the  thermal  noise  at  voltages  just  below  the  NDR 
region  increases  with  well  number  because  the  slope  of  the  I-V  curve  gets  steeper.  For 
structure  A1  the  thermal  noise  is  practically  negligible,  whereas  for  the  triple-well 
structure  A3  the  thermal  noise  actually  becomes  dominant  (Fig.  20).  Hence,  from  a 
quantitative  standpoint,  the  present  simulation  results  indicate  that  multiple-well  resonant 
tunneling  structures  may  not  have  a  better  overall  noise  performance  than  their  single¬ 
well  counterpart.  However  we  still  think  the  simulation  results  are  promising  because  the 
present  model  does  not  account  for  scattering  in  the  well(s).  Inclusion  of  scattering  would 
broaden  the  overall  transmission  function  and  thereby  smear  out  the  sharp  increase  in  the 
current  that  occurs  when  the  transmission  function  enters  the  Fermi  sea  on  the  emitter 
side.  Hence,  both  the  conductance  G  and  the  thermal  noise  contribution  to  the  noise 
spectral  density  would  decrease.  On  the  other  hand,  we  believe  the  shot  noise  reduction 
effect  would  not  be  significantly  affected  by  scattering.  By  our  model,  it  is  essentially  an 
electrostatic  effect,  dependent  primarily  on  the  amount  of  charge  stored  in  the  well(s)  and 
not  on  the  nature  of  the  tunneling  process.  This  is  supported  by  previous  experimental 
measurement  of  significant  shot  noise  reduction  in  single  and  double-well  RTDs  in  which 
sequential  tunneling  was  dominant  [32]-[34]. 

The  investigation  of  shot  noise  in  multiple-well  InGaAs/AlAs  structures  could  not 
be  extended  to  quintuple-barrier  structures  due  to  numerical  difficulties. 
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The  same  kinds  of  simulations  have  been  carried  out  for  Ino.53Gao.47As/InP 
structures.  InP  is  a  less  promising  material  as  a  potential  barrier  to  electrons  than  AlAs 
because  the  conduction  band  offset  with  lattice-matched  InGaAs  is  only  0.263  eV  (versus 
1 .2  eV  with  AlAs).  This  leads  to  smaller  peak-to-valley  current  ratios,  but  it  has  the 
advantage  of  being  much  easier  to  grow.  Although  no  experimental  data  are  available  to 
validate  the  model  for  this  material  system,  the  results  are  interesting.  A  standard  single¬ 
well  RTD  and  three  graded  quantum  well  structures  have  been  investigated.  The  design 
of  these  structures  is  shown  in  Fig.  21  and  the  simulation  results  are  presented  in  Figs.  22 
and  23.  From  Fig.  24,  the  same  decrease  of  the  minimum  shot  noise  factor  with  the 
number  of  well  is  obtained  as  in  Fig.  19.  However  it  seems  the  curve  flattens  at  about 
four  wells.  This  is  probably  due  to  the  fact  that  the  potential  barrier  is  not  very  high  for 
this  material  system,  so  that  electrons  are  able  to  transit  elastically  over  the  top  of  the  last 
barrier  of  the  biased  quadruple-well  structure.  Quantization  of  electrons  in  the  fourth  well 
is  thus  reduced,  and  the  presence  of  this  well  does  not  assist  the  shot  noise  reduction 
effect. 


structure  B1  Structure  B2 


Bias  Voltage  (V)  Bias  Voltage  (V) 

Fig.  22  Simulation  results  for  structure  B1  ((a),  (b)  and  (c))  and  structure  B2  ((d),  (e)  and  (f)) 


2.6.  Promising  Application 

From  a  practical  standpoint,  the  main  advantages  of  multiple-well  resonant 
tunneling  diodes  may  be  their  low  noise  characteristics,  as  analyzed  here,  and  their 
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broader  NDR  region  than  in  single- well  devices  [3 1],  A  rough  estimate  of  the  speed  can 
be  obtained  by  calculating  the  differential  capacitance  C  =  dQj dV  and  resistance 
K  =  dVjdl  of  the  device.  The  maximum  frequency  of  operation  around  a  given  bias  can 
then  be  estimated  by: 

/  = — — 

InRC 

For  the  triple-well  InGaAs/AlAs  structure  (structure  A3,  Fig.  18)  at  the  bias  for  which  the 
shot  noise  reduction  effect  is  maximum  (0.61  V),  computation  and  differentiation  of  the 
total  negative  charge  and  differentiation  of  Fig.  18  leads  C  «  8.0x10  *  F-cm  ^  and  R  = 
6.2x10'^  Q-cm^,  corresponding  to  an  RC  time  constant  of  0.5  ps  and  a  maximum 
frequency  (27tRC)'’  of  320  GHz.  Although  this  is  less  than  the  fmax  of  the  best  single- well 
RTDs,  it  is  still  in  the  competitive  arena  of  modem  high-speed  devices. 

The  main  application  envisioned  for  the  multiple-well  resonant  tunneling 
structure  is  the  fabrication  of  a  low  noise  microwave  vertical  field-effect  transistor  using 
the  tunneling  stracture  as  the  channel  of  the  device.  Many  groups  reported  the  fabrication 
of  resonant  tunneling  transistors  using  different  techniques  [35]-[40].  The  most  common 
design  is  to  add  a  Schottky  gate  on  the  side  of  the  resonant  tunneling  diode.  Its  principle 
of  operation  is  to  control  the  tunneling  current  by  modulation  of  the  depletion  region 
formed  around  the  gate.  The  extent  of  this  depletion  region  determines  the  effective  area 
through  which  electrons  can  tunnel.  Using  a  triple-well  structure  for  the  channel  of  the 
device  would  allow  achieving  a  better  noise  figure  without  deteriorating  its 
transconductance.  However  the  scaling  of  the  tunneling  current  with  the  effective  area  of 
the  device  is  not  perfect,  especially  at  room  temperature.  Many  challenges  in  the 
fabrication  of  resonant  tunneling  transistor  have  to  be  overcome  before  it  is  a  practical 
device. 


3.  Resonant-Piasmon-Mediated  Nanostructures 

Besides  quantum-well  heterostructures,  another  interesting  way  to  get  charge 
storage  in  semiconductors  is  embedded  metallic  or  semimetallic  layers,  or  similar 
nanoparticles.  One  of  the  more  interesting  means  of  forming  bulk  metallic  defects  in  GaAs 
has  been  low-temperature  growth  LTG  (~200°C)  whereby  excess  as  atoms  are  incorporated 
uniformly  in  the  material  during  molecular-beam  epitaxy  (MBE)  [41].  By  an  annealing 
step,  many  of  the  excess  As  atoms  precipitate,  creating  embedded  semimetallic 
nanoparticles  but  maintaining  high  crystallinity  in  the  GaAs.  In  addition.  As-related 
defects  capture  most  of  the  free  carriers,  increasing  the  room-temperature  dark  resistivity 
well  above  the  values  of  normal  undoped  GaAs. 
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'able  I.  ErAs:GaAs  Material  Characteristics 

Sample 

ErAs  layer  thickness 
(monolayer) 

Fill  Fraction  (%) 

1 

0.8 

0.56 

2 

1.2 

0.85 

3 

1.6 

1.13 
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More  recently,  the  useful  characteristics  of  the  LTG-GaAs  have  been  obtained  at 
normal  growth  temperature  (535°C)  by  incorporation  of  ErAs  during  the  growth  [42].  Like 
excess  As  in  LTG  GaAs,  the  ErAs  self-assembles  into  nanometer-scale  semimetallic 
particles  that  accumulate  excess  carriers  and  act  as  fast  recombination  centers.  The  ErAs 
approach  offers  better  control  over  the  transport  characteristics  (e.g.,  the  photocarrier 
lifetime)  of  the  material.  Fortunately,  during  the  course  of  the  present  project  the  group  of 
Prof.  Art  Gossard  at  the  University  of  California,  Santa  Barbara,  made  various  samples  of 
ErAs: GaAs  available.  In  addition,  they  prepared  samples  of  ErSb:GaSb  -  a  more  recent 
nanoparticle  materials  system  with  even  better  materials  characteristics  than  the 
ErAs:GaAs.  Samples  from  both  materials  systems  were  studied  optically  and  electrically. 

3.1  Optical  Characterization  of  ErAs:GaAs 

In  this  research  we  carried  out  the  first  known  measurements  of  the  optical 
attenuation  of  the  ErAs:  GaAs  material  at  photon  energies  well  below  the  GaAs  band  edge. 
Three  samples  were  measured  at  room  temperature,  each  consisting  of  30  layers  of  thin 
ErAs  separated  from  each  other  by  40  nm  of  undoped  GaAs.  As  listed  in  Table  I,  the 
samples  differed  only  in  the  thickness  of  the  ErAs  layers  of  0.8, 1 .2,  and  1 .6  monolayer  for 
samples  1 , 2,  and  3,  respectively.  All  three  were  grown  by  MBE  under  conditions 
identical  to  those  described  previously  except  for  the  growth  temperature  of  630°C. 

Studies  of  the  microstructure  in  samples  grown  at  535°C  show  that  the  ErAs  forms  nm- 
sized  islands  [43]  about  3  to  4  ML  (0.9  to  1.2  nm)  high  and  about  1  to  2  nm  in  lateral 
dimension.  They  form  isolated  islands  at  low  ErAs  deposition  and  network-like  structures 
at  higher  deposition.  Unpublished  studies  of  samples  grown  at  630°C  (separate  from  the 
samples  studied  here)  show  that  the  deposited  ErAs  forms  isolated  islands,  which  are  again 
about  3  to  4  ML  (0.9  to  1 .2  nm)  high,  but  now  extend  over  a  broad  range  up  to  ~1 0  nm  in 
the  lateral  dimensions  [44].  The  island  density  at  a  deposition  of  1.6  ML  is  approximately 
4x10”  cm’^.  The  changing  deposition  is  expected  to  mainly  affect  the  island  density. 

Optical  transmission  measurements  on  the  three  samples  were  earned  out  with  two 
instruments:  a  Shimadzu  Model  UV-3101PC  grating  monochromator  (PbS  detector  and 
tungsten  source,  2  cm”  resolution)  operating  between  0.8  and  2.5  micron,  and  a  Jasco  Model 
420  Fourier-transform  infrared  (FTIR)  spectrometer  (pyroelectric  detector,  incandescent 
source,  4  cm”  resolution)  operating  between  2.5  and  10  pm  wavelength.  Pieces  of  samples 
1 , 2,  and  3  were  mounted  in  metal  holders  and  the  transmission  coefficient  was  determined 
by  dividing  the  transmitted  spectrum  h  by  the  incident  spectrum  h  as  measured  by  the 
instruments  and  recorded  in  units  of  absorbance,  A  =  -logic  (Vh)-  To  separate  out 
attenuation  effects  in  the  epitaxial  layer  from  the  500-pm-thick  semi-insulating  (SI)  GaAs 
substrate  below  it,  the  absorbance  of  a  bare  SI-GaAs  substrate  from  the  same  lot  was  also 
measured. 

Because  at  all  measurement  wavelengths  there  was  no  evidence  of  standing  waves 
or  other  complicating  multiple-pass  effects,  we  applied  Beer's  law  in  conjunction  with 
single-pass  transmission  through  each  semiconductor  interface  to  write  Tsampie  =It/Ii  “ 
TiT2exp(-aepitepi)  exp(-asubtsub)  where  Ti  is  the  power  transmission  coefficient  through 
the  first  (air-epitaxial)  interface,  T2  is  the  transmission  through  the  second  (substrate-air) 
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Fig.  23.  Normalized  intensity  transmission  coefficient  vs.  wavelength  for  three  ErAs:  GaAs  samples  at 
room  temperature.  The  precipitous  drop  in  transmission  below  0.9  micron  is  caused  by  the 
fundamental  band-gap  of  GaAs.  The  weak  attenuation  features  around  1 .8, 2.4,  and  6.5  pm  are  caused 
by  H2O  vapor,  and  the  sharper  feature  at  4.3  micron  is  attributed  to  CO2 


interface,  aepi(tepi)  is  the  attenuation  coefficient  (thickness)  of  the  ErAs-bearing  epitaxial 
layer,  and  asub(tsub)  is  attenuation  coefficient  (thickness)  of  the  GaAs  substrate.  In  each 
sample  the  substrates  were  double  side  polished  so  that  the  air-GaAs  interfaces  were  of 
comparable  optical  quality  and  Ti  wTa  =  T.  Hence  the  ratio  of  the  sample  (substrate  + 
epitaxial  layer)  transmission  to  the  substrate  transmission  (i.e.,  the  normalized 
transmission)  can  be  expressed  as  Tsample  /Tsubstrate  =  [T^exp(-aepi-tepi)exp(-asub-tsub)] 
/[T^exp(-CXsub  tsub)]  ~  CXp(-(Xepi  tepi)-  5  ^Itd  thuS  CCepi  ~  -ln{  Tsample /Tsubstrate  ‘l/tepi 

The  lack  of  multiple-pass  interference  effects  in  our  samples  is  attributed  to  the  use  of 
incoherent  radiation  in  both  spectrometers,  and  to  mounting  the  samples  at  the  point  of  an 
uncollimated  beam.  Note  that  the  largest  transmission  measured  in  the  GaAs  substrate  was 
approximately  0.55  in  good  agreement  with  a  single-pass  approximation  It »  T^  and  a  T  from 
Fresnel’s  (normal  incidence)  equation  of  4ns/(ns+l)^  =  0.71  for  ns  =  3.3  (accepted  value  for 
GaAs  in  the  near-IR  at  Xs  well  beyond  the  band  edge).  This  implies  that  ttsub  is  negligible  in 
our  measurement  wavelength  range  -  a  plausible  result  since  the  present  substrates  are  semi- 
insulating  (very  low  free  carrier  absorption). 

Plotted  in  Fig.  23  is  the  experimental  normalized  transmission  measured  through 
Samples  1-3  between  X  =  0.5-  and  10  pm.  To  maximize  data  quality,  the  signal  from  0.5  to 
2.5  micron  is  taken  from  the  grating  spectrometer,  and  the  signal  from  2.5  to  10.0  micron  is 
taken  from  the  FTIR  spectrometer.  The  dominant  features  in  each  spectrum  are  the 
precipitous  GaAs  band  edge  at  -0.85  micron  and  a  broad,  asymmetric  attenuation  feature 
centered  around  2.48, 2.55,  and  2.59  micron  for  Samples  1, 2,  and  3,  respectively.  The  depth 
of  the  attenuation  feature  grows  continuously  between  Samples  1  and  3.  The  three  curves 
also  display  spurious  absorption  bands  centered  around  2.6,  and  6.5  microns  that  are 
attributed  to  residual  H2O  vapor  in  the  spectrometer,  and  a  sharp  feature  at  4.3  micron  that  is 
attributed  to  CO2. 
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Fig.  24.  Plots  of  experimental  attenuation  coefficient  vs.  wavelength  for  the  three  ErAs:  GaAs 
samples  at  room  temperature.  The  attenuation  length  is  assumed  to  be  the  entire  1 .2-pm 
epitaxial  layer  of  each  sample.  Also  shown  are  the  theoretical  curves  fit  to  experiment 
according  to  the  Maxwell-Garnett  and  Drude  models  described  in  the  text. 


Plotted  in  Fig.  24  is  the  attenuation  coefficient  derived  from  the  normalized 
transmission  using  the  expression  given  above.  Note  that  our  experiments  cannot 
distinguish  attenuation  by  absorption  from  attenuation  by  scattering.  Since  the  epitaxial 
thickness  is  1200  nm  for  each  sample,  the  peak  attenuation  coefficient  of  the  2.5-micron 
attenuation  feature  is  approximately  1300, 1850,  and  2380  cm  *  for  Samples  1,  2,  and  3, 
respectively.  The  asymmetry  of  this  feature  appears  clearly  as  a  long-w^avelength  tail  that 
produces  easily  measured  attenuation,  even  at  X,  =  10  microns.  The  peak  attenuation  values 
are  surprisingly  high  considering  that  the  material  presumably  responsible  for  the  attenuation, 
the  ErAs  layers,  constitute  only  about  1%  of  the  epitaxial  volume. 


3.2  Modeling  of  ErAs  Nanoparticles:  Surface  Plasmons 


To  understand  the  experimental  results  we  have  considered  several  possible 
attenuation  mechanisms  including  free-carrier  optical  effects,  quantum-size  effects,  and 
photoionization  from  the  ErAs  nanoparticles.  The  most  accurate  mechanism  has  been  by 
far  the  free  carrier  effects,  which  are  predicted  to  be  very  strong  because  of  a  surface 
plasmon  (Frohlich)  resonance  in  the  ErAs  nanoparticles.  Our  free-carrier  optical  model 
treats  the  ErAs:  GaAs  epilayer  as  a  composite  medium  consisting  of  independent  ErAs 
nanoparticles  separated  by  insulating  GaAs.  Consistent  with  previous  experiments  [45] 
and  theory  [46],  the  ErAs  is  treated  as  a  semimetal  using  the  Drude  model  to  account  for 
the  free-carrier  optics.  As  such,  we  can  write  for  the  dielectric  function  of  the  ErAs  [47] 


Of, 
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where  Sb  is  the  (real)  background  dielectric  constant,  (o  is  the  circular  optical  frequency, 
®P  is  the  plasma  frequency,  and  y  (=  l/x)  is  the  free-carrier  scattering  frequency  and  x  is 
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the  scattering  time.  In  the  composite  medium,  the  GaAs  matrix  is  expected  to  modify  the 
dielectric  function  in  a  manner  that  depends  critically  on  the  size  and  shape  of  the  ErAs 
nanoparticles  [48].  In  our  samples  and  over  the  experimental  wavelength  range,  we 
know  that  the  ErAs  particles  are  much  smaller  than  X.  If  we  also  assume  that  the  ErAs 
particles  occupy  a  relatively  small  fraction  of  the  total  material  volume  and  are  optically 
independent,  we  can  apply  one  of  the  truly  venerable  theories  of  optical  physics,  the 
Maxwell-Gamett  (M-G)  formulation.  In  its  general  form,  the  M-G  predicts  a  composite 
dielectric  function  of  the  form 

_  fP^S  ^2) 

f  fp 

where  Em  is  the  dielectric  constant  of  the  GaAs  matrix,  f  is  the  volumetric  fill  fraction  of 
the  ErAs  component,  and  p  is  a  numerical  factor  that  depends  on  the  shape  and  size  of  the 
ErAs  nanoparticles.  For  the  special  case  of  spherical  particles,  the  volume  dependence 
vanishes  and  p  reduces  to  3sm/(ss  +  2  em);  leading  to  the  expression. 
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From  lattice-images  taken  of  material  similar  to  that  measiued  here,  we  know  that  the 
ErAs  nanoparticles  are  closer  to  ellipsoidal  in  shape  than  spherical,  but  we  proceed  with 
the  spherical  assumption  because  of  its  great  simplicity  and  utility  in  providing  a  first 
explanation  of  the  experimental  results. 

To  apply  Eqns  (1)  and  (3)  in  the  present  samples,  we  must  know  values  for  5 
material-dependent  quantities.  The  volumetric  fraction  is  known  accurately  in  each 
sample  from  MBE  thickness  measurements  (RHEED  oseillations)  and  is  listed  in  Table  I. 
The  GaAs  matrix  has  a  dielectric  constant  that  is  slowly  varying  over  the  experimental 
range  with  an  average  value  of  em  »  10.9.  To  get  ©p,  y,  and  eb  in  the  nanoparticles,  we 
must  rely  on  theoretical  estimation  since  the  semimetallic  nature  of  ErAs  means  that 
electrons  and  holes  will  both  contribute  to  any  transport  effect  (including  plasma 
excitations),  and  contribution  of  each  carrier  type  to  the  plasma  frequeney  will  depend  on 
the  detailed  band  structure.  It  is  known  that  ErAs  has  a  rocksalt  band  structure  so  that  the 
electrons  reside  in  ellipsoidal  valleys  characterized  by  a  longitudinal  and  transverse 
effective  masses,  mt  and  mj,  respectively.  The  holes  reside  in  three  warped-sphere 
bands,  the  light-hole,  heavy-hole,  and  split-off  bands  having  effective  masses  mLH>  hihh, 
and  mso,  respectively.  Fortunately  a  great  deal  of  experimental  work  based  on  the 
Shubnikov-de-Haas  effect  coupled  with  band  structure  calculations  has  also  produced  the 
concentrations  of  each  carrier  type,  as  listed  in  Table  II.  Given  these  quantities  and  given 
the  applicability  of  the  semiclassical  model  of  transport  of  independent  carriers,^  an 
interesting  exercise  in  deductive  reasoning  is  to  show  that  the  plasma  frequency  of  a 
many-carrier  solid  is  given  by  [49]  (in  MKS  units), 

(O  = - r - = - (  — 

So  ^0  W 


37 


where  e  is  the  electron  charge,  eq  is  the  vacuum  permittivity,  N  is  the  number  of  relevant 
bands,  and  the  symbol  <p/m>  symbolizes  the  averaging  process  defined  by  (4)  over  the 
densities  and  effective  masses.  For  example,  in  the  four-band  model  of  ErAs  we  find: 


2  ^ 
=- 

€, 


Pe^lh^hh^so  +  Plh^e^hh^so  +  PhhVlh^so  +  Pso^e^lh^hh 


(5) 


Given  the  material  characteristics  in  Table  II  this  results  in  a  linear  plasma  frequency 
(cop/27i)  of  6.7x1 Hz,  corresponding  to  a  plasma  (cutoff)  wavelength  of  0.45  pm.  This 
is  considerably  longer  than  the  plasma  cut-off  in  most  metals  since  the  total  carrier 
concentration  in  the  ErAs  (8.3x10^°  cm‘^)  is  so  much  lower  (e.g.,  in  copper  p  =  8.4x10^^ 
cm'^). 

The  final  two  quantities  in  the  model  are  the  least  certain,  but  are  also  not  as 
crucial  to  our  interpretation  of  the  experimental  data.  The  scattering  frequency  is  a 
complicated  function  of  particle  size  and  shape,  but  simplifies  to  y  =  ybuik  +  3<VF>/2d  for 
spherical  particles,  where  ybuik  is  the  bulk  scattering  frequency,  <vf>  is  the  average  Fermi 
velocity,  and  d  is  the  diameter  of  the  sphere  [50] .  ybuik  can  be  estimated  from 
measurements  made  on  the  conductivity  of  continuous  ErAs  films  on  GaAs  at  room 
temperature,  a 1.5x10'*  /Q-cm  [51].  The  bulk  scattering  can  then  be  estimated  through 
the  following  variation  of  (4), 


y  = 


(6) 


where  y  is  assumed  to  be  independent  of  carrier  type.  Substitution  of  the  experimental 
conductivity  along  with  the  value  e^<p/m>  =  1.58x10^*^,  yields  ybuik  ~  1.05x10*'*  s  *,  or 
Xbuik  =  9.9  fs  -  a  reasonable  result  since  it  is  so  close  to  the  Drudian  scattering  time  in 
common  metals  (e.g.,  x  «  10  fs  in  Al).  To  get  the  surface-scattering  contribution,  we 
write  the  Fermi  velocity  in  each  band  according  to  the  well-known  expression 

V/, ,  =  (^/ffj,.)(3;^  ^nd  calculate  an  average  by<  >=  ’• 

This  leads  to  the  result  for  ErAs  of  <vf>  =  1 .14x10^  m/s  -  a  value  comparable  to  that  in 
common  metals  (e.g.,  vf  in  Al  =  2.02x10^  m/s),  in  spite  of  the  much  lower  concentration 
of  free  carriers  in  ErAs.  For  the  ErAs  nanoparticles  previously  imaged  by  TEM  methods 
the  size  of  the  particles  is  in  the  range  of  1  to  2  nm.  When  combined  with  the  above 
value  of  <vf>,  the  estimated  range  of  y  is  1.7x10*^  to  8.5x10*'*  s'*  (x  =  5.8x10  *^  s  to 
1 .2x1  O'*  ^  s),  well  below  the  bulk  value.  Because  we  do  not  yet  know  the  exact 
distribution  of  sizes  and  shapes  of  the  nanoparticles  in  the  present  samples,  the  diameter 
will  be  treated  as  a  fitting  parameter  and  compared  with  the  above  range  for  credibility. 


3.3  Surface  plasmon  effect  on  optical  transmission 


The  final  physical  quantity  in  the  model  is  eb  -  the  background  dielectric  constant 
associated  with  the  core  electrons.  As  in  many  metals  this  quantity  is  hard  to  measure 
and  calculate,  so  we  will  assume  it  to  be  1.0.  In  the  wavelength  range  of  measurement, 
the  dependence  on  Eb  is  weak  anyv^'ay  because  both  the  real  and  imaginary  part  of  the  e  is 
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determined  by  the  large  GaAs  background  sm  and  plasma  effects.  Given  the  5  quantities 
and  the  dielectric  function  of  (3),  we  can  calculate  the  attenuation  and  propagation 
coefficient,  a  and  P,  respectively,  from  classical  electromagnetics, 

a,p  =  (2(o/c){(p'  8V2){[1+(s"/£’)^]’^^  -,+  1 } } (7) 

where  c  is  the  speed  of  light  in  vacuum  s’ (s”)  are  the  real  and  imaginary  parts  of  the  M- 
G  dielectric  function  of  (3),  and  p’  is  the  relative  magnetic  permeability  [52].  Although 
ErAs  is  known  to  have  a  large  positive  p’  owing  to  unfilled  inner-shell  paramagnetism, 
such  effects  are  almost  certainly  relaxed  at  the  present  wavelengths,  so  that  we  can 
assume  p’=l.  Note  that  we  call  this  an  absorption  rather  than  attenuation  coefficient 
since  the  M-G  theory  treats  the  entire  epitaxial  layer  as  a  uniform  composite  medium. 
Backscatter  is  also  accounted  for  by  this  formalism  but  shows  up  through  modification  of 
p  (i.e.,  the  real  part  of  the  refractive  index)  and  a  modification  of  the  reflection  coefficient 
via  the  Fresnel  equations. 

Fig.  24  shows  the  computed  value  of  the  intensity  a.ttenuation  constant  2a 
according  to  the  quantities  calculated  for  each  of  the  three  ErAs:GaAs  samples.  The 
only  fitting  parameter  -  the  diameter  of  the  spherical  particles,  d  -  is  set  equal  to  1 .2  nm. 
Remarkably,  the  theoretical  curves  display  very  similar  behavior  to  the  experimental 
ones.  Each  theoretical  curve  shows  an  asymmetric  resonance  peak  with  a  long- 
wavelength  tail  that  has  significant  absorption  even  at  X-  =  10  pm.  The  absorption 
increases  monotonically  in  peak  strength  with  increasing  ErAs  fill  fraction,  although  the 
theoretical  peak  wavelength  for  each  sample  -  Xp  =  2.1  pm  -  is  somewhat  shorter  than  the 
experimental  peaks.  This  disagreement  in  peak  wavelength  causes  a  significant 
discrepancy  between  experimental  attenuation  and  theoretical  absorption  on  the  short- 
wavelength  side.  Note,  however,  that  the  slope  of  the  theoretical  and  experimental  curves 
is  nearly  the  same.  Note  also  that  the  fit  value  of  d  =  1 .2  nm  is  in  the  range  of  expected 
nanoparticle  thickness  (0.9  to  1 .2  nm)  in  the  present  samples,  but  is  well  below  the 
expected  lateral  dimension. 

To  further  support  the  absorptive  model  we  also  computed  the  M-G  intensity 
reflection  coefficient  R  according  to  (3)  and  (7)  and  the  Fresnel  equation,  R  =  1(  - 

l)/(s’^^  +l)f  •  We  find  no  resonance  in  R  around  2.5  pm,  although  there  is  a  weak 
variation  vs  X  consistent  with  the  Drude  model.  To  explain  the  discrepancies  between 
model  and  experiment,  we  note  that  a  more  refined  analysis  would  represent  the 
nanoparticles  not  as  perfect  spheres  but  as  ellipsoids  or,  possibly,  spheroids.  In  either 
case,  the  M-G  and  full-wave  scattering  theories  would  both  predict  the  resonant 
absorption  to  shift  to  longer  wavelengths  relative  to  the  peak  of  the  spherical-particle 
curve  [48].  Because  an  optical  sum  rule  applies  to  the  free-carrier  optics,  the  peak 
absorption  would  also  diminish  relative  to  that  of  the  spherical-particle  model,  and  thus 
broaden  the  overall  signature  -  more  consistent  with  experimental  results. 

Given  the  qualitative  success  of  the  model,  it  is  compelling  to  elucidate  the  nano¬ 
optics  behind  the  effect.  It  is  well  known  that  metal  particles  dispersed  at  low  fill 
fraction  in  an  insulating  host  can  display  resonant  attenuation  below  the  bulk  plasma 
frequency  associated  with  normal  electromagnetic  modes  of  the  particle.  When  the 
normal  mode  is  quantized,  these  become  the  so-called  surface  plasmons.  Remarkably, 
the  fundamental  mode  can  be  predicted  from  the  M-G  theory  in  the  limit  of  very  small  fill 
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fraction.  In  this  case  (3)  displays  a  pole  in  the  dielectric  function  at  the  frequency  where 
Ss «  -2£m  a  condition  credited  to  Frohlich  [53].  For  the  Drude  model  of  free-carrier 
optics,  this  condition  becomes  simple  to  interpret  in  the  range  y  «  ©.  Then  from  (1),  the 
frequency  that  satisfies  ss  «  -2sm  is  simply  to  =  (Op/(8B  +  28m)*^=  top  [54],  often  called 
the  Frohlich  resonant  frequency.  In  the  present  samples  with  cop  =  4.2x10^^  Hz,  A,p  = 

0.45  pm,  8m  =  10.9  and  8b  «  1.0,  we  find  top  =  8.8x10''’  Hz  and  X.p  =  2.1  pm,  in  good 
agreement  with  the  peak  of  the  model  attenuation  curves. 

Note,  however,  that  in  the  present  samples  we  do  not  satisfy  y  <<  to.  Rather, 
because  of  the  small  size  of  the  nanoparticles  and  the  associated  surface  scattering,  our 
estimation  of  y  yields  a  value  of  9.6x1 0''’  s''  -  slightly  larger  than  top.  Hence,  the  simple 
interpretation  given  above  is  not  strictly  valid.  Nevertheless,  the  model  resonance 
persists  and  its  peak  remains  near  ^p  =  2.1  pm.  Although  it  is  not  simple  to  prove  this 
mathematically  from  (1),  (3),  and  (7),  it  can  be  argued  from  the  following  classical 
viewpoint.  The  normal  mode  that  gives  rise  to  the  surface  plasmon  and  Frohlich 
resonance  is  a  result  of  classical  electromagnetics  and  the  Drude  model.  In  other  words, 
it  is  just  one  form  of  a  Lorentzian  atomic  resonator  with  zero  resonance  frequency.  As 
such,  when  the  damping  factor  (y  in  our  case)  increases,  the  resonance  broadens  and 
weakens  but  in  a  gradual  fashion,  so  the  peak  persists  and  the  peak  frequency  remains 
approximately  constant. 

It  is  important  to  contrast  our  results  with  the  previous  work  by  Nolte  et  al.  who 
sought  to  observe  Frohlich  resonances  in  metal  precipitates  in  GaAs  and  AlGaAs 
samples,  particularly  As  in  low-temperature-grown  GaAs  [55,56].  First,  it  is  possible 
that  their  measurements,  which  were  only  taken  out  to  1 .5  pm,  did  not  go  far  enough  in 
wavelength  to  observe  a  resonance  in  the  As;  GaAs  samples.  And  the  presence  of 
standing  waves  in  the  AlGaAs  layers  of  the  As:  AlGaAs  samples  would,  in  our 
experience,  make  it  very  difficult  to  observe  the  resonances  at  the  attenuation  strength 
observed  in  our  samples.  Second,  it  is  quite  possible  that  our  ErAs  nanoparticles  are 
much  more  monodisperse  in  size  and  shape  than  the  As  precipitates  in  the  low- 
temperature-grown  samples.  More  lattice-imaging  experiments  will  be  required  to  prove 
this.  But  if  true,  we  anticipate  that  because  of  the  strength  of  the  surface-plasmon 
resonance,  ErAs  nanoparticles  will  offer  a  new  means  of  engineering  the  optical 
properties  of  GaAs  and  InGaAs  at  near-infrared  wavelengths  where  widespread 
technological  development  continues  for  optical-fiber  telecommunications  and  other 
applications. 


Table  II.  Band  Properties  of  ErAs^ 

Band 

Carrier 

Concentration  (cm'^) 

Effective 

Mass  (xme) 

Electron 

4.14E+20 

0.21 

Light  Hole 

1.19E+20 

0.28 

Heavy  Hole 

2.93E+20 

0.46 

Split-Off 

2.00E+18 

0.15 
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Conclusion 


Through  this  research,  a  formalism  that  enables  the  study  of  the  noise  properties 
of  multiple-well  resonant  tunneling  diodes  has  been  developed.  The  shot  noise  reduction 
effect  is  modeled  as  the  result  of  the  electrostatic  repulsion  encountered  by  the  injected 
carriers  due  to  the  charge  stored  in  the  quantum  well(s).  Simulation  results  are  compared 
successfully  to  experimental  results,  and  then  carried  out  for  the  first  time  on  multiple- 
well  resonant  tunneling  structures.  It  appears  that  multiple-well  structures  show  superior 
noise  properties,  i.e.,  a  greater  shot  noise  reduction  effect,  than  their  single-well 
counterpart.  These  simulations  results  must  be  confirmed  experimentally  but  are  already 
promising  for  the  use  of  multiple-  quantum-well  structures  in  new  low-noise  devices  and 
circuits. 

In  the  spirit  of  “Fermi-electronics”,  this  research  has  also  spawned  the  discovery 
of  very  strong  plasmon-mediated  optical  absorption  in  a  semimetal-semiconductor 
composite;  ErAs:GaAs.  Strong  absorption  signatures  have  been  observed  just  below  the 
GaAs  band-edge,  and  have  been  explained  theoretically  in  terms  of  a  surface  plasmon 
resonance.  This  research  paves  the  way  for  a  class  of  optically-  and  (perhaps)  electrically 
triggered  switches  in  which  the  gating  function  is  carried  out  by  resonant  plasmons. 
Although  similar  in  spirit  to  plasma-electronic  devices  recently  studied  by  M.  Shur  et  al., 
resonant  plasmon  coupling  is  much  stronger  so  that  the  devices  can  have  much  higher 
electrical  or  optoelectronic  gain. 
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